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ABSTRACT 

We calculate the higher order moments in a sequence of models where the initial den- 
sity fluctuations are drawn from a xt distribution with a power-law power spectrum. 
For large values of v the distribution is approximately gaussian, and we reproduce the 
values known from perturbation theory. As v is lowered the distribution becomes pro- 
gressively more non-gaussian, approximating models with rare, high-amplitude peaks. 
The limit v = 1 is a realization of recently proposed isocurvature models for producing 
early galaxy formation where the density perturbations are quadratic in a gaussian 
field. 
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1 INTRODUCTION 
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The standard paradigm for the formation of large scale 
structure is that quantum fluctuations during an inflation- 
ary epoch seeded initially small, gaussian fluctuations in 
density, which grew through the action of gravitational in- 
stability in a universe whose dominant constituent is cold £jv 
dark matter (CDM). This theory has proved very predic- 
tive and agrees with a wide range of observational data, as 
however several of the assumptions are difficult to test with 
high precision. The assumption that the initial fluctuations 
have a gaussian distribution holds more generally than just 
in the inflationary CDM model, as the central limit theo- 
rem implies that fluctuations emerging from many uncorre- 
cted, random processes will be nearly gaussian. Though it 
is indeed a plausible assumption, and a prediction in the 
simples t inflationary models, it has received only limited (IPeebles 1980| : 



gaussian random field has been studied extensively (see 
Strauss & Willick ( [l995|) for a review). Of particular interest 
are the Sn, defined in terms of the volume averaged corre- 
lation functions (cumulants of the probability distribution 
function) 



d n d 3 r 2 ■ ■ ■ d 3 r N £(ri, 



(1) 



(2) 



In the mildly non-linear regime the growth of gaussian initial 
conditions by gravitation al instabil i ty predicts that t he Sn 
are independent of scale (Fry 1984; Bernardeau 1992) for a 
scale-free spectrum. For a top-hat filter of the density field 
and a pure power-law spectrum, P ( k) oc fc" , it can be shown 
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Heavens 199£) and large-scale structure 


(Bouchet ct al. 1993 
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1994 
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Feldman, Kaiser & Peacock 1994; Stirling & Peacock 


1996 


Colley 1997). 


In the latter case the situation is made 



Colombi 1993; 



Goroff et al. 1986 



Juszkiewicz, Bouchet & 



B^rn*ard*eluTT994^ that to lowest non-trivial 



order in perturbation theory (£2 *C 1), 



more difficult by the action of gravity which turns an ini- 
tially gaussian random field into a non-gaussian field once 
the modes become non-linear. 

Fo r gaussian fluctuations the only non-trivial moment 
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with a very weak de pendence on the matter density, f2o 
(Bouchet et al. 1995). These results have been confirmed 
by N - body simulations ( Coles fc Frenk 1991 1_ Bouchet et al 



form tl ej30vrerj3p_cctnm^^£/j^l4_^^ ct al. 1994; Bernardeau 



1 the higher-order moments of the field carry addi- 1994 ; [Lokas ct al. 1995 ; [Baugh, Gaztanaga fe Efsatathiou 



gaussia 

tional iiformation. The evolution of the higher order mo- 
ments induced by gravitational instability in an initially 
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1995; |Juszkiewicz et al. 1995; Colombi et al. 1996[) a nd are 



is the 2 point correlation function, or its Fourier Trans- 1992 1_ Weinberg fc Cole 1992j; Lahav et al. 1993 ; |Juszkiewicz 



known to be insensitive to redshift space distortions ( Lahav 
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et al. 1993; |Colombi et al. 1996| ; |fflvon et al. 1995|) on the 
mildly non- linear scales where we will be working (see ^). 

Measurements of the Sn have been performed on 
several galaxy catalogues in two and three dimensions; 
see e.g. Table 1 in Hui & Gaztanaga (199£) or Kim & 
Strauss ( jl998| ). We expect that the Anglo- Australian Two 
Degree Field survey (2dFQ) and the Sloan Digital Sky Sur- 
vey (SDSS0) will provide excellent catalogues for estimation 
of the Sn for different sub-samples to high order. Assum- 
ing growth through gravitational instability from initially 
small gaussian fluctuations one can use the additional infor- 
mation contained in the Sn to relate the properties of the 
observed galaxy distribution to the fluctuations in the un- 
derlying dark matter. The observations of CMB anisotropies 
have enhanced our faith in gravity as the engine of growth. 
The extra assumption required in this step is that the initial 
conditions were gaussian. 

As with other assumptions in the standard paradigm, 
our assumption of gaussianity for the initial fluctuations 
should be tested against observations. One difficulty with 
testing for non-gaussianity has been the lack of a simple, 
predictive theory which describes what form, out of the infi- 
nite possibilities, the non-gaussianity should take. While in- 
flationary models generically predict gaussian fluctuations, 
not all inflationary models are "generic" in this sense. And of 
course there are other, non-inflationary, models of structure 
formation which predict non-gaussian perturbations. The 
best known of these non-gaussian theories are models based 
on topological defects. However they have problems fitting 
the observational data and are difficult to model properly 
requiring expensive simulations even for the initial condi- 
tions. 

Moments have been predicted for a variety of non- 
gaussian models under various approximations. Results for 
the unsmoothed moments 5*3 and S4 for arbitrary initial 
conditions in perturb a tion theory have been present ed by 
( |Fry fc Scherrer 1994 jChodorowski fc Bouchet 1996] ). Ap 



pro ximations to defect models have been studied by (Jaffc 



Gaztanaga & Mahonen 1996 ; Gaztanaga & Fosalba 



1994 

199S ) . The latter authors also discussed an approximate 
calculation of the low-order moments for the model of Pee- 
1 999b|) , which shall discuss further below. Wein- 
1992) studied models obtained from non-linear 



bles ( |1999a 
berg & Cole 

mappings of initially Gaussian fields. Per haps the closest an 
cestor to this work, however, is that of (Coles et al. 1993) 



who used N-body simulations to calculate moments for a va- 
riety of non-gaussian models, including a xi model similar 
to that developed below. 

In this paper we investigate the predictions for the low- 
est moments, S3 and S4, in a simple, parameterized, non- 
gaussian "model" . The model has useful interpolating prop- 
erties. It has one parameter, v. The limit v — > 00 recovers 
the gaussian result, with lower values of v being progres- 
sively more non-gaussian. Low values of v mimic cosmolog- 
ical models with rare, high-amplitude peaks. 

The outline of the paper is as follows: in the next section 
we discuss the non-gaussian model in more detail, in §[] we 
discuss our technique for calculating the moments S3 and 54 



* http:/ /meteor. anu.edu.au/~colless/2dF 
t http://www.sdss.org 



using an N-body code, plus the tests we have run. In §^ we 
summarize the results, and in §rawe discuss the implications. 



2 THE NON-GAUSSIAN MODEL 



The COBE team (Kogut et al. 1996) first introduced a non- 
gaussian model where the real and imaginary parts of the 
Fourier- (or in their case 1-) space perturbations, are drawn 
from independent xi distributions with v degrees of free- 
dom, adjusted to have zero mean and scaled to give the right 
power spectrum. While the precise physical significance of 
this model is unclear, it serves as a well defined reference 
with useful interpolating properties, as the xi distri bution 
becomes m ore gaussian as v — > 00. It was shown in (Kogut 
et al. 1996) that the COBE data prefer v — > 00, with the 



gaussian model being 5 times more likely than any other 
model tested. 

We have modified this model to specify xi distributions 
in real rather than fc-space. We proceeded as follows: Let 
(j>i(x) be independent, gaussian random fields of zero mean. 
Define, for integer v, 



(5) 



Then ip will be xi distributed, in real-space, with zero mean. 
For v — > 00 the fluctuations will be gaussian. As v is lowered 
the model becomes progressively more non-gaussian, with 
more "rare" peaks. The limiting case v — 1 has the per- 
turbations quadratic in a gaussian field, much as in recently 
prop o sed non-gaussian inflation models ( Linde & Mukhanov 



1997 



Antoniadis, Mazur fc Mottola 1997| ; |Peebles 1999a 



Peebles 1999b). The reduced moments of the initial condi- 
tions are easily calculated for a xi distribution. The mean 
is v, and 

((X 2 ~v) 2 ) = 2v, (6) 
((X 2 -^) 3 ) = 81/, (7) 
<( x 2 -^) 4 ) = 48^ + 3(2^) 2 . (8) 



A final feature of this model, not shared by other pro- 
posals in the literature which involve non-linear mappings 
of Gaussian density fields, is that the distribution of fluctu- 
ations will remain approximately the same when smoothed 
on a variety of length scales. This is related to the scaling 
of this model: {ip n } scales as {^ 2 ^™^ 2 . We show this feature 
explicitly for the case v = 5 in Fig. where the distribution 
of initial fluctuations from one of our simulations is given for 
different smoothing scales, spanning a factor of 5 in scale. 

It is straightforward to generate initial conditions for 
this model with a power-law spectrum. For i = 1, • ■ ■ , v we 
generate independent realizations of a gaussian field 4>i(x), 
with correlation function £i(r), on a grid in the usual way. 
We then square <f>i and sum over i to obtain tp, which we 
use as an initial field to generate displacements using the 
Zel'dovich approximation as described below. 

The only remaining step is to choose the correlation 
functions £i(r) so that the final correlation function £$(r) 
has a given form. Writing 6$ = $ — (ffi) and using the 
fact that 4>i is a gaussian one has 



(9) 
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where A 2 
fields d>i 



fc 3 P^(fc)/(27T 2 ) 

This gives £<j>{r) = 



4 



Figure 1. The distribution of the initial density fluctuations for 
v = 5, smoothed with a gaussian filter of R = 1, 3 and 5 grid 
cells in a 64 3 simulation. Note that the distribution is skew on all 
scales. 




k (h/Mpc) 



Figure 2. The power spectrum of the initial conditions for v = 1 
generated by the method described in the text. The solid line is 
the n = —1 target spectrum and the squares indicate the average 
power spectrum from 50 realizations of the ICs. The deviation is 
less than 1% over more than an order of magnitude in scale. 



Note that since the <f>i are independent, the correlation func- 
tions for the components of ij) simply add. For simplicity we 
choose them to be equal. 

It is important to notice that £^ > for all separations, 
so the "integral constraint" (i.e. the integral of f (r) over all 
space must vanish) cannot hold for this spectrum. However 
we shall be interested in power-law spectra for which £ > 
anyway, so this technical point will not affect us. We shall 
choose the power law index of the 4> field = —2 so that 
— 2ri0 +3 = —1, close to the slope observed for CDM 
models on the scales of interest. In the absence of finite box 
size corrections, we get = — 1 by choosing 



4W - 1 (f 



is the power spectrum of the 
> /2)A 2 ,{k f r)- x or 



(11) 



We set k, = 0.2/iMpc" 1 and adjust A* so that A^(fc*) = 
0.25, or as — \. The same normalization is used in the test 
simulations described below. 

We have checked whether the resulting power spectra 
for our non-gaussian models match the = — 1 form we 
expect. Since the evolution is a non-linear function of the 
initial conditions, it is important to do this comparison on 
the initial conditions rather than the evolved spectrum. We 
find that there is a departure from a pure power-law at both 
high-fc and low-fc. This is to be expected: A 2 is given by an 
auto-correlation of AJ, and in the simulations A^, is miss- 
ing modes with wavelengths longer than the box or smaller 
than the mesh scale. We have corrected for this effect by 
modifying the input spectrum of the <f> field to ensure an 
initial power-law spectrum for A 2 . While A 2 changes by 
more than 3 orders of magnitude over the range of scales we 
simulate, the correction factor we apply always lies between 
1 and 3 (in power), with the largest correction at low-fc. 
Note that we are not directly modifying the initial power 
spectrum of the ip field here - it is the convolution integral 
in the <f) field that requires this numerical correction. 

To sum up: The procedure outlined above provides an 
algorithm for generating realizations of a random field whose 
1-point function is a yt (as has been explicitly checked in 
the simulations, see §f|) and whose 2-point function has a 
power-law spectrum to within 1% over more than an order 
of magnitude in scale. We show this explicitly for v = 1 
in Fig. ^. This initial realization of the tp field can then 
be evolved using N-body simulations to study the effect of 
gravitational growth on these initial conditions, as described 
in the next section. 



3 CALCULATING THE MOMENTS 

We use iV-body simulations to calculate the moments 53 and 
5*4 in the quasi-linear regime. A discussion of several of the 



relev a nt numerical issues can be found in ( luszkiewicz et al 



1995 ; |Colombi et al. 1996: |Baugh, Gaztanaga fc Efsatathiou 



1995|; Pzapudi fc Colombi 1996j ; |Kirn fc Strauss 1998| ; flam 



& Bertschinger 1995). We use a particle-mesh (PM) code, 



AJ(*) = A2(£) 



(10) 



described in detail in Meiksin, White & Peacock (1999). 
We have simulated critical density universes (fio = 1) with 
power-law spectra in boxes of size 150/i _1 Mpc < Lbox < 
250/i" 1 Mpc on a side, so that the fundamental mode was 
always well in the linear regime and the box represented a 
"fair" sample of the universe. To minimize finite volume ef- 
fects, we will work on scales less than about 0.15i?box, where 
^box = (4vr/3)i? box (see below). 

We have run simulations with either 64 3 or 128 3 par- 
ticles and a 64 3 or 128 3 force "mesh". All the simulations 
were started at 1 + z = 20 to 30 and run to the present 
(z = 0). The evolution was done in log of the scale fac- 
tor a = (1 + z)~ . The time step was dynamically chosen 
as a small fraction of the inverse square root of the max- 
imum acceleration, with an upper limit of Aa/a = 4 per 
cent per step. This resulted in a final particle position er- 
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ror of less than .1 pe r cent of the box size. As described 
in Meiksin et al. ( 199E ) , the code reproduces the non-linear 
power spectrum excellently as compared wi th se mi-analytic 
fitting functions (e.g. Peacock & Dodds (1996)) or other 
N-body codes. The initial conditions were generated from 
an initial density field, with either an A^art or 7V mGS h FFT, 
using the Zel'dovich approximation. The particles were ini- 
tially placed either on a uniform grid, or at random within 

,1/3 



m Peacock & Dodds ( 199q ). Tests indicate the most reli- 
able conditions are when the resolution imposed by finite 
particle number matches the force resolution from the mesh 
and when the particles were displaced from random posi- 
tions within cells, i.e. the co nditio ns used in the simulations 
described in Meiksin et al. (199E). 



At the end of each simulation 8, 192 sphere centers were 
thrown down at random within the volume and the number 
of particles in each of multiple concentric top-hat spheres 
computed. The radii of the cells was constrained to be larger 
than ~ 2.5 (mesh) cells and smaller than 0.2i?box- Since the 
smallest sphere is more than 5 grid cells across, the effects 
due to the finite resolution of the PM code are minimal. An 
explicit check with the higher resolution simulations verified 
this expectation, though to achieve per cent level accuracy 
we need to work on scales larger than 4 mesh cells (8 in 
diameter) for the larger boxes where the grid scale is less 
non-linear. This is also to be expected: in the non-linear 
regime the structure on the grid scale is determined mainly 
by the collapse of larger wavelength modes, which should be 
well ev olved b y the PM code at all times ( Bouchet, Adam 
& Pelht 1985). Thus as the grid scale becomes more non- 
linear, the effects of the finite force grid on perturbations a 
few grid cells across are lessened. We have been careful to 
ensure that the non-linear scale at z = has a wavelength of 
several mesh zones. Comparison with different box sizes and 
resolutions suggested the larger radius cells, with radii ap- 
proaching 0.2i?box, were affected by the finite box size, with 
Sa being more affected than S3 . For this reason we restricted 
the largest cell in the analysis to be less th an 0.15-Rbox- An- 
alytic arguments (Hui & Gaztanaga 1998) suggest that our 



boxes are large enough to avoid finite volume biases at the 
few per cent level on these scales. 

The moments of the counts-in-cells distribution were 
calculated and averaged over many (> 100) simulations with 
different realizations of the same initial spectrum and statis- 
tics. Writing 



N- N 
N 



(12) 



where N = (TV) and (■ • ■) represents an ensemble average 
(which w e approximat e by an average over the simulations) 
we have ( Peebles 198(i| ) 



H2 
IM 



TV 

S^N' 1 + 2/V" 2 



3/4 + nfj, 2 N- 2 



67V" 3 



(13) 
(14) 
(15) 



Our estimates of £jv are used to calculate 53 and S4, with the 
errors (including correlations) estimated from the run-to-run 
scatter and propagated in the usual manner. Note we aver- 
age the \iu over the simulations and then compute the Sn 
rather than averaging the Sn computed from each simula- 




20 40 60 

Simulation Number 

Figure 3. Estimates of S^(R = lO/i" 1 Mpc) from the cells 
thrown for each of a set of n = —1/2 simulations as described 
in the text. The average of these values is the dashed line, and 
the S4 computed from the entire suite of simulations before av- 
eraging is the rightmost (solid) point with error bar. This value 
agrees well with the perturbation theory expectation, shown as 
the solid line. 



tion. This avoids biases introduced because (x/y) 7^ (x)/{y) 

From the work of Hui & Gaz- 



( Hui fc Gaztanaga 199 ^ 
tanaga ( |1998] ) we estimate that their so called "estimation 
biases" are 0(1%) for the cases of interest here. We show in 
Fig. |i] the estimates of S4 from the cells thrown for a set of 
simulations with gaussian initial conditions and power-law 
index n — —1/2. The average of these values is shown as 
the dashed line and the 54 computed from the entire suite 
of simulations as described above is the final (solid) point 
with error bar. This value agrees well with the perturbation 
theory expectation - the solid line - whereas the naive aver- 
age does not. In this Figure the error bars are overestimated 
for each realization by the small number of cells thrown, but 
this error is reduced by the large number of realization used. 

An alternative to the counts-in-cells method described 
above, is to perform a maximum lik elihood fit to the entir e 
distribution function of the counts (Kim & Strauss 1998). 



This method minimizes errors due to finite volume effects 
and shot-noise. However it requires one to know a priori 
a functional form for the distribution func tion. The Edge- 
worth expansion, used by Kim & Strauss (1998), is a valid 
expansion only near the peak of the distribution. To re- 
move unph ysica l oscillations it must be regularized. Kim 



& Strauss ( 1998 ) remove the oscillations by convolving the 
Edgeworth expansion with a Poisson distribution, to account 
for the very sparse sampling they performed. However once 
the number of particles in a cell becomes appreciably greater 
than unity this convolution no longer regularizes the distri- 
bution and a different expansion is needed. Typically models 
for the distribution of counts involve many higher order mo- 
ments, beyond S4, making implementation of this method 
difficult. For these reasons we have chosen to use the more 
traditional moments method. Another alternative to throw- 



ing cells (Lokas et al. 1995) is to interpolate the final density 
field on a grid, and compute moments of the smoothed field 
by summing over the grid points (the smoothing is done us- 
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Figure 4. The simulation results for S3 and S4 for 3 power-law 
spectra with n = —1 (solid squares), n = —1/2 (open circles) 
and n = (open triangles) as a function of £2- The perturba- 
tion theory results, valid for £2 <C 1, are shown as the horizontal 
dashed lines. Error bars (shown for S4) are calculated from the 
run-to-run scatter in the simulations. For 53 the la error bars are 
approximately the same size as the points themselves. The error 
bars on the points are correlated. 



ing FFT methods). For our implementation - using top-hat 
cells and small grids - this method, while significantly faster, 
was not as accurate. For fine enough grids and smooth in- 
terpolation (e.g. CIC or TSC) we would expect to obtain 
the same results as the direct counts-in-cells method, but 
we did n ot inv estigate this in detail. Finally the method of 
Szapudi ( 199cj ) for throwing an effectively infinite number 
of cells was too CPU intensive. For the low order moments 
we are considering, the traditional method allows sufficient 
cells to be thrown that Szapudi's method is not required. 

As a test we first simulated gaussian initial conditions 
with pure power-law spectra with n = — 1, n = —0.5 and n = 
0, for which the results are known analytically (see Eqs. || 
^) when I2 <C 1. The agreement between the numerical and 
analytic results, shown in Fig. ^, was very good, showing 
that numerical effects were under control. It is interesting 
to note that for n ~ the perturbation theory results are 
obtained only when £2 is quite small, in contrast to the case 
where n ~ — 1, for which the perturbation theory results 
are a good approxi mation even for £2 ^ 0.1. This agrees 
with the results of ( Colombi et al. 1996 ), see their Fig. 5. 
For n = — 1 (our fiducial model) the numerical results and 
analytic predictions agreed to better than 5 per cent, for 
both S3 and S4, when £2 <§C 1. This was the case for either 
gaussian initial conditions, or runs with the non-gaussian 
model of with v 1. This is an impo rtant result, since i t 
argues that the simulations reported in (Meiksin et al. 199£ ) 
can be safely used to estimate correlations bet ween power 
spectru m b ins in the mildly non-linear regime (Meiksin & 
White " 



99£ ) 



Having established that the code performs as expected, 
we simulated models with = — 1 as described in §^ for 
varying values of v. The results are described in the next sec- 
tion. Before leaving these tests however, we show in Fig. ^ 
the values of 53 and S4 for the n = — 1 spectra with gaus- 



Figure 5. The value of 53 and 54 from the n = —1 simulations 
in real space for a critical density universe (solid squares), an 
fim = 1 — f2/Y = 0-3 universe (open circles) and in redshift space 
in the A model (open triangles). 



sian initial conditions, but varying Q m and in redshift rather 
than real space. Our redshift space results use the distant 
observer approximation, i.e. for each of the simulations we 
added the z-component of the velocity, in units of the Hubble 
constant, to the z-component of the position before throw- 
ing the cells. These results confirm that redshift space dis- 
tortions and variations in £l m have a negli gible impact upon 
the values of ^3 and 1 S4 that we calcula te ( Lahav et al. 1993 ; 
Hivon et al. 19951). 



Colombi et al. 1996 



4 RESULTS 

Our main results are shown in Fig. ^. Here we plot 
S 3 (10/i _1 Mpc) and S4(T0 h^Mpc) as a function of v. We 
have chosen this scale because, for our normalization, £2 — 
0.1 at lO/i Mpc. For n = —1 and gaussian initial condi- 
tions we would expect S3 ~ 2.857 and 5*4 ~ 13.89 in the 
weakly non-linear regime. These values are shown as the 
dotted lines in Fig. ^. Note that 53 and S4 increase rapidly 
as the model is made more non-gaussian. 

Once the model is non-gaussian we no longer have rea- 
son to expect that the Sn will be independent of scale. How- 
ever non-linear effects process the initial non-gaussianity in 
a non-trivial way, as seen in Fig. ^[ Simple scaling arguments 
from the initial conditions, see Eqs. (|l3|-|l5|), would suggest 



that £ 3 n 



and CI 



(12/u)£i. This is to be 



compared to the gaussian prediction that £| 



and 



£f auss oc Q. Defining S' N such that S' N would be indepen- 
dent of v and £2 for a xt field we have 

S' 3 = S3&v) 1/2 (16) 
S' 4 = S 4 (£» (17) 

for unsmoothed fields. For a pure xt field, S' 3 = y/E ~ 2.82 
and S4 = 12. By calculating the moments of the initial den- 
sity field from which the Zel'dovich displacements are gen- 
erated, we have explicitly checked that these values are in 
fact attained in the initial conditions on the grid (i.e. with- 
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Figure 6. The moments, S3 and 54, as a function of v for the 
non-gaussian model described in ^. These values are evaluated at 
10 h~ x Mpc where £2 — 0.1. The limit v = 1 has the perturbations 
quadratic in a gaussian field, as in recently proposed isocurvaturc 
models. The limit v — > 00 (the rightmost point) recovers the 
gaussian result. The predictions of perturbation theory assuming 
gaussian initial conditions, S3 ~ 2.857 and £4 ~ 13.89, are shown 
as the horizontal dashed lines. 




Figure 7. The modified moments, S' 3 and S' 4 , defined in the text, 
as a function of v as in Fig. ft The dashed lines show the values 



2.82 and Si 



12 of the x?, initial conditions. Note that 



J 3 — ^■"^ "4 — x ^ ul ullc "X-v 
gravitational instability has induced a dependence on v. 



out any smoothing) . Any deviation is at the 1 per cent level 
for S'3 and S 4 over a wide range of v. We show the evolved 
values of S3 and S 4 vs v in Fig. (?], with the gaussian point at 
v — 00 omitted. Notice that gravity has modified the initial 
conditions so that the S' N are no longer independent of v. 
The curves rise to larger u, indicating that the gravitation- 
ally induced skewness and kurtosis are a larger fraction of 
the initial conditions as the model becomes more gaussian. 

In addition, the effects of non-linearity change the scal- 
ing with £2- For example, for the v = 1 model the Sn are 
much less dependent on £2 than the S' N , varying by 0(10%) 
for 0.03 < £2 < 0.3. This behaviour can be explained as the 



gravitational contribution becoming more important than 
the initial condit ions a s clustering evolves, as anticipated by 
Fry & Scherrer ( 1994 ). This result suggests that one must 
use care in inferring from the scale-independence of the Sjv 
that the initial fluctuations are gaussian. To properly make 
this inference it may be necessary to use very large scale or 
high redshift measurements. 



5 CONCLUSIONS 

We have used N-body simulations to calculate the low order 
moments of the mass density field for a sequence of non- 
gaussian models with power-law spectra. We obtain good 
agreement with the predictions of perturbation theory for 
gaussian initial conditions. As the real-space perturbations 
are made progressively more non-gaussian the moments de- 
part significantly from the values of Eqs. (|[^|) as shown in 
Fig. ^| Moreover, gravitational evolution modifies the scaling 
of the Sjv with £2 and v. For example the scale-dependence 
of the Sjv does not follow that of the initial conditions over 
the range the simulations probe: 0.03 < £ 2 < 0.3. 
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